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The density-wave structures are studied in quasi-one-dimensional atomic gas mixture of one- 
component bosons and two-component fermions using the mean-field approximation. Owing to 
the Peierls instability in the quasi-one-dimensional fermion system, the ground state of the system 
shows the fermion density wave and the periodic Bose-Einstein condensation induced by the boson- 
fermion interatomic interaction. For the two-component fermions, two density waves appear in each 
component, and the phase difference between them distinguishes two types of ground states, the 
in-phase and the out-phase density- waves. In this paper, a self-consistent method in the mean-field 
approximation is presented to treat the density-wave states in boson-fermion mixture with two- 
component fermions. From the analysis of the effective potential and the interaction energies, the 
density-waves are shown to appear in the ground state, which are in-phase or out-phase depending 
on the strength of the inter-fermion interaction. It is also shown that the periodic Bose-Einstein 
condensate coexists with the in-phase density-wave of fermions, but, in the case of the out-phase 
one, only the uniform condensate appears. The phase diagram of the system is given for the effective 
coupling constants. 

PACS numbers: 03.75Kk, 03.75.Ss, 71.45Lr, 75.30Fv 



I. INTRODUCTION 

Since the experimental achievement of Bose-Einstein condensates (BEC) 0,0 and Fermi-degenerate 
states |3| , the researches on the ultra-cold atomic gas have much stimulated the quantum many-body 
physics; it can create an artificial many-body system with the controllable the atom number, trapping 
potential geometry, the interatomic interactions et.al. 

In recent years, the quantum-degenerate atomic gas of the boson-fermion mixture has been realized 
with the advanced cooling technique (the sympathetic coolin g ) | 4l |5 | , and a lot of experimental and 
theoretical studies have been done for it [rjR 11 H ITol HI ITaTmTl4|. 

The quasi-one-dimensional atomic (Q1D) gas of the boson-fermion mixture should open up many 
interesting possibilities in the quantum many-body system 0]: the Tonk- Girardeau gas 0,0], and 
other regimes of trapped ID bosons 0] and of fermions 0]. As discussed in 0|, the Q1D atomic gas 
can be realized, for example, in the axially-deformed harmonic oscillator (HO) trap with the axial and 
radial frequencies Lo a , r , which satisfy Lu a <C uj r and should be larger than the interatomic interaction 
energies and the healing length of the system (TEl I20L l2l| . Experimentally, the quasi-one-dimensional 
systems has already been realized in the boson atomi c gas la . |23. 12a . l24j , where the BEC solitons 
propagating in the axial direction have been observed |25|, I26U27I]. 

A particular property of the fermionic Q1D system is the occurrence of the density wave due to the 
Peierls instability with the frequency of twice the Fermi wave number 2fcp. In |15| . we investigated 
the ground-state properties and the collective excitations in the mixture of one-component boson and 
one-component fermion, and showed that the fermion density wave appears in the ground state of the 
system, which is similar with the charge density wave (CDW) in the Q1D electron system pi I2H30I ] . 
and we also found that the BEC of the bosons becomes spatially periodic (periodic BEC). 

In the present paper, we investigate the ground state properties of the Q1D atomic gas at zero- 
temperature (T = 0) of one-component boson and two-component fermion, which interact through 
the atomic collisions. In this system, two density waves can appear in each component of fermions in 
the ground state, and two types of the ground-state structure are distinguished with the relative phase 
between them, the in-phase and the out-phase density- waves (Fig. 1); the latter corresponds to the 
spin density wave (SDW) in the electron system |3C| . The density- wave structure in the ground state 
depends on the strength of the interatomic interactions; especially important is the inter-component 
interaction of fermions, which is absent in the previous study |15|. Also, the structure of BEC, uniform 
or periodic, should be influenced by the fermion density-wave structure. 
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FIG. 1: In-phase and out-phase density waves in fermion densities Pf(z) (solid curve) and pi(z) (dashed 
curve) . 

In Sec. II, we introduce a model Hamiltonian for the system, and propose a new method to treat 
the periodic order parameters for boson and fermion-density fields, which are proved to be self- 
consistent in the mean-field approximation, and construct the effective potential with this method. In 
Sec. IH, the density-wave structures in the ground state are classified by analytical evaluation of the 
effective potential for the cases of the repulsive/attractive fermion- fermion interactions; the results are 
summarized in the phase diagram for the effective coupling constants of the interatomic interactions. 
The final section is devoted to the summary of the paper. 



We consider the Q1D system of atomic gas mixture of one-component boson and two-component 
fermion at T = 0. In this paper, we call two components of fermion "up" and "down", and denote 
them with the symbol f and J,. It should be noted that they should not only be the up and down 
components of spin 1/2 atom, but any two of hyperfine levels or different atomic species. 

The masses of the bosons and the fermions are denoted by nib and m/. In this paper, we consider the 
case of degenerate two-component fermion for simplicity, and assume that the up and down fermions 
have the same mass (also, we assume that the other properties, the interaction strengths et.al., are 
degenerate) . 

The boson-fermion mixture is assumed to be trapped in the HO potentials with the same frequencies: 



If the HO potentials are enough long in the z-direction, u> r 3> u> n , the system becomes Q1D. (For the 
other conditions in the realization of the Q1D system, see jl5L |20L l21| . ^ 

In the Q1D system, the dynamical degrees of freedom of atoms are represented by the ID field 
operators, 4>{z) for boson and ipt t ±(z) for fermions. Using these field, the Hamiltonian of the system 
is given by 



II. FORMULATION 



A. 



Hamiltonian 



U bJ = m bJ {cu 2 (x 2 + y 2 ) + u? a z 2 } /2, 



(1) 



H 




' h 2 d 2 

2mb dz 2 





+ U' f (z)-n f Mz) 





(2) 
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where fj,bj are one-dimensional chemical potentials for the bosons and the fermions. The coupling 
constants gbb,bf,ff m eq. (J2J) represent the strengths of the boson-boson, boson-fermion and fermion- 
fcrmion interactions, respectively. In the case of the large attractive boson-boson interaction (gbb <C 0), 
the BEC instability appears and competes with the Peierls instability. In the present paper, to 
concentrate on the density wave states, we consider the case of the repulsive boson-boson interaction 
(gbb > 0) only. For the Q1D atomic gas, the coupling constants are determined by the (3D) s- 
wave scattering lengths abb.bf.ff of the atomic collisions |3l|: gbb — 2TiLo r abb, gbf — 2hui r abf, and 
gff = 2huj r a,ff, where uj r is the radial frequency of the (3D) HO potential The U' h j(z) in 

eq. (|21) are the z-directional part in the original HO potential QJ, which are neglected in the present 
calculation. (About the effects of the HO potential, see 0|.) 



B. Density- wave ground state 



It is well known in the Q1D fermion system 0, |2^, [2t| [3(J and also in the boson-fermion mixture 
flU that the homogeneous state is unstable (Peierls instability) for the formation of the density wave 
with the wave number 2k F . To describe the density- wave states, we introduce the 2&F-periodic order 
parameters for the boson and fermion fields: 



(<f>(z)) = $(z) = b Q + 2b c cos(2k F z), 

Pl (r) = {tp\(z)i> l (z)} = f + 2f c cos(2k F z + ^)+2f b cos{2k F z), 
p T (r) = {4>\{z)^{z)) = fo + 2f c coa(2k F z + <p + 9) + 2f b co8(2k F z), 



(3) 
(4) 
(5) 



where (p is the phase-difference angles between boson- and fermion-density waves, and 9 is that between 
the up- and down-fermion density waves; 9 = (n) corresponds to the in-phase (out-phase) density 
wave of fermions (Fig. 1). The b c and f c in eq. (JSJ are the amplitudes of the boson and fermion 
density waves, and the third components with the amplitude fb in have to be added to reflect 
the fermion-density waves induced by the periodic BEC (the term with the amplitude b c in the boson 
order parameter). Thus, the amplitude fb will be shown to be proportional to the b c at the end 
of this section. The constants bo and fo represent the homogeneous terms in the boson and fermion 
densities, which are determined by the boson and fermion average densities, rib and rif after the spatial 
integration of the local densities: 



n b = b\ + 2b 2 c , n f = n/, T = n fd = f 



(0) 



In the present study, we do not consider the case where the spin-polarized states with n/,-f — n/,j 7^ 
become stable; in the cases of the very large fermion density (or the very strong fermion-fcrmion 
interaction), such a polarized state might be stable (the collective ferromagnetic states) |3^ . 
Using eqs. (I3I5II for the Hamiltonian J2J, we obtain the mean- field Hamiltonian density: 



it-MF = T~(-F + rii 



(7) 



where the fermion part TL F includes the contributions from the fermion kinetic term and the fcrmion- 
fermion and boson-fermion interaction terms, and the boson part He comes from the boson kinetic 
and boson-boson interaction terms. Using the periodic regularization with the spatial length L and 
taking the limit L — > 00, the fermion part for the momentum —2k F < k < 0, denoted by Tip, becomes 



H F = 



dk 
2^ 

dk 

9ff ,2 1 
2 9ff 



iP ] (k + 2k F ) 
^(fc + 2fc F ) 



t r 



A'" 



e(k) 



e(k) 



A b + A' b 
e(A:4 

A fc 4 



-2k F 

-a; 

e(k 



^ ^{k + 2k F ) 

2k F ) 



^i(k + 2k F ) 



A) cos(6») 



2A/A|,cos ( - 



^2 



Af 



(8) 



The fermion part, 7i F , corresponding to the positive momentum < k < 2k Fl is obtained by the 
replacement k F — > ~k F in eq. JSJ), and the total fermion part becomes 7i F = 7i F + TCp. The e(k) is 



the single-particle energy e(k) = H k /2rrif — fif +3///0 + gbfnb, and the gap functions A are defined 

by 



A b = 2b b c g bf , A f = f c g ff , A' b = f b g ff . 



(9) 
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FIG. 2: Illustration of the quasi-particle spectrum, E\. ^ (k). We scaled quantities by the fermi energy and 
fermi momentum, and set A|,|/2eF = 0.1. 



Using the similar technique with the Bogoliubov transformation, the Tip is diagonalized (for the 
detailed calculation, see Appendix A): 



dk 

2^ 



9ff f 2 _ 1 

2 h 



J2 [E-(k)a-\k)a-(k)+E+(k)at\k)at(k) 



9ff 



A) cos(6>) + 2A f A' b cos (^j 



cos <p 



(10) 



where a^(fc) are the annihilation operators for the quasi-particle with the momentum —2kp < k < 
in the energies: 



Ef{k) = [e(fc) + v F {k + k F )] ± ^v 2 F {k + k F f + |A T | 2 , 
Ef(k) = [e(fc) + v F (k + k F )} ± ^(fc + fc F ) 2 + jA^ 2 , 



(11) 
(12) 



where vf is the Fermi velocity. The quasi-particle energies for the momentum < k < 2kF are ob- 
tained by replacing kF, Vp — > — kp, —vp in eqs. (|11I12|) . The effective gap functions Af ^ in eqs. (|11I12|1 
are defined by 



IAtI 2 



A' b + e^A 



/I 



(A 6 



,'\2 



2(A b + A' b )A/ cos(<^) + Aj, 



|A X | 



A b + A' b + e^+^Af ' = (A b + A' b f + 2(A b + A ft )A/ cos(0 + <p) + A* 



(13) 
(14) 



In Fig. |2 we show the quasi-particle energy in eqs. (|11I12|) . It should be noted that the energy gaps 
appear at k = ±kp in the quasi-particle spectrum. 

In the density-wave ground state, the up- or down-fermions occupy the lower branch of the quasi- 
particle spectrum so that the energy density of fermions in this state is given by the summation of 
E^^(k) in eqs. (|11I12I) up to the Fermi momentum fc F : 



(H F ) 



dk 
2^ 



[E~(k, -k F )9(k F - k)9{k) + E~(k, k F )0{k F + k)9(-k)] 



-9fff ~ — 
9ff 



Aj cos(6») + 2A/A' fc cos ( - I cos ( <ys + - 



A? 



+ 9fffo - M/ + 9bfn b 



V± V 

2 ^ 

»=T4 L 

-9ff ft ~ ^~ 
9ff 



A 2 2e F + v/44 + Af 



A 2 , cos(6>) + 2A/A' 6 cos - cos Us + - 



A b 2 



(15) 
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where we have used (H F ) = (7~t F ) + (T~L F ) = 2 {Hp). 

Now we turn to the evaluation of the mean-field energy for the boson part TLb in eq. JSJ) . Substituting 
eq. into the boson part, we obtain the boson energy density up to 0(b 2 ): 

(H B ) = 2[e b (2k F ) + 6 9 fbl]bl + 9 fbt + 6 9 fbi 



e b (2k F )+5^b 2 



,2 , .966 ,2 



2e b (2k F )b 2 + 10 



9bb &b , 9bb u2 



2W f ' ~ b ° nb 



9bb 2 



kp 
mbUb 



9bb 



2 

Ah x 2 

9b f , 



(16) 



In the derivation of the last line, we have used the solution of n b = b 2 l + 2b 2 (the average boson density) 
and A & = 2b b c g bf : 



bl = - 



26; = - 



ft* > 2 



5b/ 



'« 2 -2[ — 



A, 

9bf 



21* 

,9bf 



(17) 

(18) 
(19) 



where we have assumed the small modulation in the boson density b 2 3> 2b 2 , and taken the solution 
satisfying 6q = rib an d 2b\ = at A& = 0. 

Thus the total energy density is given by the sum of eqs. (|1 5I1C3I) : 



£ tot {A b ,A f ,A' b ;<p,9) = (H MF ) = {H F ) + (H B ). 



(20) 



Explicit representation of the energy difference £tot(A b , A/, A' b ) between the density-wave state and 
the normal state (A;, = A/ = A/, = 0) is given by 



£ tot (A b ,A f ,A' b ) = £ tot (A b ,A f , A' b ) - £ tot (Q, 0, 0) 



+2nf6 F 



A a 17 A ? , 2e F + y/4e F + A 2 S 



( k 



\m b n b 



9ff 



A? cosfl + 2A/A^cos ( - ] cos ( <p + - ) + A' 2 



l-21og 



cos6< + 2 



A^ 

A/A^ 
(2 £F )2 



A^ 
2eir 



2 

1 / Afc_ 



cos - COS </? + - 



.966 



2e F 



At 

56/ 



(21) 



(22) 



where dimensionless parameters £ and £ (effective coupling constants) have been introduced: 



_ 9bf 



Trhg bb v F 

= 9ff 
Trhv F ' 



/ v F m 



\v B m b 



(23) 
(24) 



and vb = \f <?66^6 / Tn b is the bosonic velocity. 

The stable ground state of the system is determined from the minimum point of eq. (|22[) . 



G 



C. Consistency of the order parameters in mean- field approximation 

Here we check the consistency of the Fermion density ansatz (0J|5j in the mean-field approximation, 
where the fermion densities p^i(z) are represented by 

Ps {z) ee (ipt(z)Mz)) 



dk 

-(^(fc)^ s (fc)) 



— (^ s (k + 2k F )Mk)y 2kFZ +C.C, (s=tl)- (25) 

In the quasi-particle representation (shown in Appendix A), the expectation values appeared in eq. (|25[1 
become 

dk 

-(^(fc)^(fc)) = n f , (26) 

^t(k + 2k F )Mk)) = 2 f g — " As (27) 

J - kF Zn 2 v /{e(fc) - e(k + 2k F )} 2 + \A S \* 

Using the above equations in eq. (|25J) . we obtain the quasi-particle representation of the fermion 
densities: 

p^( z )= n f + 2 I ^ / " lutT " i; n cos(2fc F z) 

:Cos(2fc F z + <^), (28) 





dk 




—(Aft + A' b ) 




/-fc F 


2lT \ 




-e(k + 2k F )} 2 ^ 


"|A T | 2 


f 


dk 




-A/ 






2lT \ 


/{<k) 


-e(fc + 2A: F )} 2 -t 


-|A T | 2 


f 


dk 




-(A 6 + A 6 ) 




)-k F 


2lT \ 


/{<*>) 


-e(k + 2k F )} 2 ^ 


-|A|| 2 


f 


dk 




-A/ 




l-kw 


2tt 


l{<k) 


-e{k + 2k F )} 2 4 


-|AJ 2 



Pl(z) = n f + 2 I — — = cos(2fc F z) 

:COs(2A: F z + + <£>). (29) 

Here, the densities /Jf,j.(^) above show the same z-dependence with those in eqs. JHEJ: it proved the 
consistency of the present method. 

In the present paper, we consider the case of n/j = i.e., no polarization, and have assumed 
the same amplitudes in ansatz (|4I5[1 . so that we can put |AjJ = |A||. From eqs. l|13lHl)l for the gap 
functions Af^, we obtain the relation between the phase-difference angles: cos(ip) = cos(6 + <p). 

Furthermore, we can derive the relation among the gap functions at the stationary points of the 
mean- field energy H2U|) . From the gap equations derived in Appendix B, we obtain 

A' b = -J;A b -A / cos(£)cos(V+£) . (30) 

Using the definition of the gap functions in , it can be written as the relation of the density- wave 
amplitudes: 

h = -/ccos (31) 

where the effective boson- fermion interaction parameter G b f is defined by 

/ k 2 

2 K F 



G bf =gi f [—^+g bb ) . (32) 
1 \m b n b I 
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In the present paper, we consider the case of gbb > 0, so that the parameter Gbf is positive (Gbf > 0). 
The relation i|30|) shows that A' b is related to a linear combination of and A/ in the ground state. 
In the next section, we will show that is proportional only to the boson amplitude Ah after fixing 
the phase-difference angles in the ground state. 



D. Interaction energy and phase differences 



To determine the ground state, we should fix the phase-difference angles ip and 9 so as to gain the 
largest interaction energy for the fixed gap functions. For this purpose, we reevaluate the ip- and 
^-dependent terms in l|22|) . included in the boson- fermion and fermion-fermion interaction terms in 
eq. J2J), in the mean- field approximation. 

Using the parametrization for order parameters (|3I4I5I) . the boson- fermion interaction term becomes 

9bf / dz ( f ) H z ) ( t ) ( z )^t( z )'>Ps(z) 
L dz 



9bfJ Q T £* t (*)<K*X^(*W.(*)>li 



2g b f 
2g b jn b nf 



n b fo + 4&o& c / c cos l - ] cos ( tp- 



46 6 c / fc 



A? 
2 b 

G b f 



(33) 
(34) 



where, in the derivation of the last line, we have used the relation (|3f Jl . which is valid for the ground 
state. Eq. I|34|) shows that the boson-fermion interaction term is independent of the phase-difference 
angles in the ground state, and plays no role in the phase difference determination. 
In the same way, we reevaluate the fermion-fermion interaction term: 



9ff / dz ^(z^iz^liz^^z) 



9ff 



dz 



2 9ff 



fo , 
2 



fc cos | - | cos I s 



+ fl 



2 9ff 



— cos - 
2 



cos [ if + - 



AG 2 



-ft 



sin ip 



(35) 



The phase-diffcrcncc angles ip and 9 should be determined to make the energy gain largest in eq. (|35|l 
for fixed values of the coupling constant Therein, the (9, </>)-dependent part appears only in the 
second term, so that the phase-difference angles are determined as the minimum point of eq. I|35l) : 
(9, ip) — (0,7r/2) for gff < 0, and (9, ip) = (tt,tt/2) for gff > 0. It should be noted that the results 
are consistent with the condition cos(ip) = cos(# -I- p) obtained in previous subsection. 



III. CLASSIFICATION OF THE DENSITY- WAVE STRUCTURES IN THE GROUND 

STATE 

Now, we analyze the minimum of the mean-field energy l|22|l with the phase difference conditions 
obtained in the last part of sec. II, to determine the stable state of the system, and show the complete 
classification of the ground state for the dimensionless coupling constants C, and £. 
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FIG. 3: Contour plots of Etot in As- A/ plane for the repulsive Fermion-Fermion interaction gff > 0. Three 
plots are for £ > £ (a), £ = £ (b) and £ < £ (c). 



A. Repulsive fermion-fermion interaction g// > 



In the case of g/f > 0, using (8, <p) = (ir, tt/2) for the phase difference angles in the ground state in 
eq. ((30(1 . we obtain A' b = — A;, (= — ^-A^. Substituting them into eq. 1(22(1 . we obtain the total 
energy density up to the order of 0(A 2 ): 



2nf€F 



1 - log 



A 2 / 




(36) 



where Abe ( 1 — 4? 



2 j A b /{2e F ), and A, = A f /(2e F ). 

In Fig.0 the contour plots of the energy density 1)36(1 are shown in Ab- A/ plane for typical values 
of the coupling constants C and £. 

Since the first term in eq. ((36(1 is symmetric with respect to the order parameters Abj, the position 
of the minimum point of eq. 1(36(1 is determined by the last two terms; the results are summarized as 



• o < e < c 



(Fig.HO, 



(X 



2C-C 



e 2 c-e ,0 



In this case, the up- and down-fermions form in-phase density wave: 



2e F ( 



cos(2kpz) 



(37) 



and the bosons make a periodic BEC: 



n b (z) ~ n b 



9bfn b 



■ cos(2fc F z) 



(38) 



which is analogous to charge-density waves (CDW) of electrons in low dimensional conductors 
[30l l33| due to the Peierls instability. We here mention the roles which the sign of g b / plays 
in the density waves. Since the amplitude of the periodic BEC in (|38|l is always positive by 
definition, the phase difference between boson and fermion density waves is determined by the 
sign of the amplitude of the fermion density wave in 1)37(1 . which is proportional to the sign of 
g b f. attractive (repulsive) boson-fermion interaction for the in-phase (out of phase). 



A} = 4e~ 



. C = £ (Fig.Eb), 

The above equation shows the existence of the rotation-like "hidden" symmetry in the ground 
state, which mixes between the bosonic and fermionic particle-hole pair order parameters. It 
also suggests the occurrence of the spontaneous symmetry breaking in the ground state and the 
existence of the Nambu-Goldstone mode as a low-lying excitation mode. 



9 



£=0.5, §=-0.6 




FIG. 4: Contour plots of £ to t in As-A/ plane for the attractive fermion-fermion interaction gjj < 0. 

• C < C < 2C (Fig. Eft) ande<2C, (A 6i A/) = (o, 2e^) . 
In this case, the fermion-fermion interaction induces the out-phase fermion- density waves: 



1 ± sin(2fc F z) 



(39) 



which forms the out-phase density wave in the two-component fermion system, analogous to the 

sdw M,M. 



When 2£ < £, we need a special care to treat the boson density wave, because the coefficient of 
A 2 term in l|36[) becomes negative, which can be traced back to the potential energy of the gap 
functions. Assembling the boson-fermion and fermion-fermion interaction terms l|34l I35fl and 
the boson energy l)lfi[) . the potential energy up to the constant, H pot , becomes 



Hpot — 



G 



bf 



1 - 



9ff 



2G 



bf 



A? 



9ff 



A 2 



(40) 



which should be negative definite for any value of the gap functions to gain the condensation 
energy. Thus, it gives A^ = for 2Gb / < gff (or equivalently for 2( < £). 



B. Attractive fermion-fermion interaction gff < 



In the case of gff < 0, using the phase difference angles (0,<p) = (0,7r/2) the relation A' b 
is obtained from eq. (|30|) in the ground state. The total energy density (|22[) becomes 



St, 



2nf€F 




An + —A 2 

2C+ICI B l£l 



/ 



(41) 



e ac+iei ] and fermion 



Fig. 0] shows a contour plot of the energy density in A^-A/ plane when ( — 0.5 and £ = —0.6. 

As is clear from eq. (|4Tj) . the coefficient of A 2 B is larger than that of A 2 , 2/(2£ + |£|) > 2/|£|, so 
that the minimum point of the energy density is always on the A b axis. 

Thus, for gff < 0, the order parameters take values (^Af,, A./j — ^2( 

densities become in-phase, p\{z) = p[(z) = rif 1 — cos(2kpz) 

those in eq. I|37|) in the repulsive fermion-fermion case gff > 0. The boson condensate has also the 
same form with (I38|l . However, the ^-dependence of the order parameter A;, is different in both cases 
as shown in Fig. [S] where the gap functions are plotted for the effective coupling constant £ when 

c = i. 



2C+ICI 

which have the same form with 
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FIG. 5: The gap functions (extreme values) at (= 1 as a function of £. 
C. Phase diagram 

Finally, in Fig. |HJ the phase diagram of the density-wave structure in the ground state is shown in 
the plane of the effective coupling constants ( and £ in eqs. 1)2^124(1 . which correspond to the gh/gbb 
and 5//. 

In Fig. there exist two phases denoted by in-phase and out-phase, which correspond to the 
fermionic in-phase density wave (9 = 0) and out-phase density wave (9 = tt) with/without the periodic 
BEC in the ground state. 

In the case of the repulsive fermion-fermion interaction (£ > 0), there exists a competition between 
in- and out-phase density waves, and the border of these two phases is given by £ = £ (rigid line in 
Fig. EJ). Another critical border line £ — 2( exists (dotted line in Fig. EJ, and, above this line, the 
ground state has the out-phase fermion density waves but no periodic BEC. 

As for the attractive fermion-fermion interaction (£ < 0), the ground state has the periodic BEC 
which is coherent with the in-phase fermion density waves. The fermion-fermion interaction affects 
the boson magnitude Ab implicitly. 

The behaviors of the order parameters and the effective gap functions I |13I14|I in Fig.Elare reflected 
in the phase diagram; A/ and A& are continuous and smooth functions of £ as stationary points of 
the energy density, and disappear at £ = and £ = 2£ respectively, as mentioned above. The line for 
the effective gap function A^j (solid line), the net energy gap at the Fermi surfaces, bends at £ = £; 
it implies the transition from the in-phase density wave to the out-phase one. 

IV. SUMMARY AND DISCUSSION 

We studied the ground state of the Q1D boson-fermion mixture of the atomic gas with two- 
component fermion and one-component boson, and showed in the consistently formulated mean-field 
approximation that the states with the periodic BEC and the fermion density waves with the wave 
number 2fcp become stable therein at T = We also showed that there exists competition between 
the in-phase and out-phase fermion density waves in the ground state, which is originated from the 
relative strength between the boson-fermion and fermion-fermion effective interactions. We clarified 
the wave-density structures of the ground state and gave the phase diagram of it. 

Finally, we make a rough estimation for the density-wave amplitudes of the atomic gas mixtures 
in the axially deformed trap potential. As possible candidates, we take rubidium isotope mixtures, 
84 Rb- 87 Rb and 86 Rb- 87 Rb, with the parameters uj a = 2tt x 10 Hz and uj r = 2tt x 15 kHz for the trap 
potential, and Nf ~ 10 3 and Nf, ~ 2 x 10 4 for the atom numbers, with which the Peierls instability 
is shown to exist in ^{|. For the scattering lengths in atomic collisions, we take the values in [3£j 
(atomic unit): 




1) 84 Rb- 87 Rb: (o//, o 6 /, a bb ) = (142, 117, 90) (C, £) = (0.397, 0.467), 



(42) 
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FIG. 6: The phase diagram of the density-wave structure of the ground state in plane. 



2) 86 Rb- 87 Rb: (a//, a bf , a bb ) = (7, 336, 90) =► (£ £) = (3.276, 0.023). (43) 

Using these parameters, we obtain the results that, in 84 Rb- 87 Rb mixture (the case 1)), the out-phase 
density wave should occur with the amplitude Af/(e F £) ~ 0.118, and, in 86 Rb- 87 Rb mixture (the case 
2)), the ground state has the in-phase one with the amplitude Af/(2e F () ~ 0.451 and periodic BEC 
with amplitudes At,/(gbfn b ) ~ 0.321. The changes of the atom number, the trap potential size and 
also the interaction strengths should make the shift of these results, and make the phase transitions 
shown in Fig. 6 expected to be observable. 



APPENDIX A: DIAGONALIZATION OF THE MEAN-FIELD HAMILTONIAN 

In this appendix, we briefly sketch the diagonalization of the mean-field Hamiltonian (jSJ under the 
periodic BEC and the fermion density waves given in eqs. © and JHJ, and show the quasi-particle 
wave functions. 

The up- and down-fermion contributions in eq. JSJ) have the same form in structure and the calcu- 
lations go completely in parallel, so that we omit the subscript of the fermion field operator in this 
appendix. The first and second terms in eq. JHJ should be diagonalized by a unitary matrix U: 



ip(k + 2k F ) 



uu- 



a (k) 
a+(k) 



e(k) A 
A* e(k + 2k F ) 



1 r E-(k) 






E+(k) 



uu- L 

a'(k) 
a+(k) 



ip(k + 2k F ) 



(Al) 



where £' ± (fc) are the quasi-particle spectrum given in eqs. (|llfl and (|12[) . and a ± (fc) is the field operator 
corresponding to them. 

The explicit form of the unitary matrix is given by 



U = 



v(k) u(k) 
u(k) v(k) 



u- 1 = U* = 



v*(k) u*{k) 
u*\k) v*(k) 



(A2) 



where 



e(k) - e(k + 2k F ) 



^/{e(k)-e(k + 2k F )} 2 + \A\< 



1/2 



(A3) 
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A /l 
AIV2 



1 - 



i{k) - e{k + 2k F ) 



{e(k)-e(k + 2k F )Y + \A\ 2 
1/2 

e(fc) - e(Jfc + 2fej?) 



1/2 



{e(fc)-e(fc + 2fc F )} 2 + |A| 2 

e(k) -e(k + 2k F ) 
{e(k)-e(k + 2k F )} 2 + \A\ 2 



1/2 



(A4) 



(A5) 



(A6) 



for — k F < k < 0. The quasi-particle wave functions, u(k) and v(k), have been normalized to satisfy 
the condition: \u(k)\ 2 + \v(k)\ 2 = \u{k)\ 2 + \v(k)\ 2 = 1. 



APPENDIX B: GAP EQUATIONS 

The gap equations of the system are given by 

I A/ cos(#) + A; cos f cos + |j } 
(Aj + Ajj cos(</?) + A/ 



4 
5// 



dk 
2ir 



2A b 



= 2 



yJ{e{k)-e(k + 2k F )} 2 + |A T 

A/ cos(y) + Ab + A b 
A /{e(fc)~e(fc + 2fc F )} 2 + |A T 

-lj A ' b+ A /C os(£)co S ^ + £ 



dk 
2^ 



dk 
2^ 



A f cos(y) + A b + Aj, 
{e(fc)-e(fc + 2fc F )} 2 + |A T 



(A b + A b )cos(y + 6>) + A / 
v /{ e (fc)- £ (fc + 2M} 2 + |A i | 2 

A/ cos(y? + g) + A b + A 6 
^{ e (fc)- e (fc + 2fc F )} 2 + |A i |^ 



A/ cos(y + 6) + A b + A' b 
{e(k)-e(k + 2k F )} 2 + \A l f 



(Bl) 



(B2) 



(B3) 



which can be derived from the energy density eq. (|21|l by taking the variation with respect to the gap 
functions Af, Ab, and A' b each other. 
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